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Abstract 

In the construction of modern Celestial Reference Frames (CRFs) the overall rotational alignment 
is only weakly constrained by the data. Therefore, common practice has been to apply a 3-dimensional 
No-Net-Rotation (NNR) constraint in order to align an under-construction frame to the ICRF. We 
present evidence that correlations amongst source position parameters must be accounted for in order 
to properly align a CRF at the 5-10 /ras level of uncertainty found in current work. Failure to do so 
creates errors at the 10-40 /ras level. 


1. Introduction to No-net-rotation Constraint Equations 

Since the adoption of the ICRF [4] by the IAU in 1997 [1], the axes which define Right Ascension 
and declination have been realized by an ensemble of quasar positions. Subsequent frames, most 
notably the ICRF2 [5], have been aligned to the original ICRF by imposing a 3-dimensional No- 
Net-Rotation (NNR) constraint in the form of the cross-product (see full derivation in appendix), 

N src 

^2 So i X As* = 0 (1) 

i= 1 

where for a set of N src sources, 

soi = (cos n 0 j cos <f 0 j, sina 0 jcosc>oj, sind'oi) (2) 

is the unit vector in the reference source direction for the i th source of coordinates (aop^o*) and 
A Si is the difference of the reference position from the estimated position. 

soi x Asi = [— cos sin d'oi cos <5oi Aaj + sinctoi A5j], (3) 

[— sin sin 5oi cos Aoj — cosao* A5j], 

[cos 2 5oi A on] 

The above triplet of constraints affects both the estimated set of on and Si and their associated 
covariance including the correlations amongst parameters. 

The MODEST linear least squares fit software [7] requires the coefficient with respect to A a* 

and ASi of each rotation constraint, Cj = (JT s o* x A si)j where j goes from 1 to 3. The result is 


dCn/dAcn = 

— cos no i sin Soi cos S a, 

dCx/dASi = 

sin aoi 

(4) 

dC 2 /dA ai = 

— sin aoi sin Sa cos <5oi , 

dC 2 /dASi = 

— COS CKoi 

(5) 

dCs/dAcn = 

cos 2 Sou 

dC 3 /dA5i = 

0 

(6) 
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Note that the constraint Eq. (1) is unweighted in the position parameters, s. The MODEST 
software implements constraints by treating them as pseudo-observations with given uncertainties. 
Specifying a vanishingly small uncertainty is equivalent to using an absolute constraint. 

Using the preceding partials (Eqs. 4, 5, 6), each Cj was separately constrained to be zero to 
within a tight constraint uncertainty of Oj = 10” 10 radians. 

2. Sensitivity of S/X, K, and X/Ka Frame Alignment to Correlations 

Our procedure was first to estimate source positions under the influence of unweighted NNR 
constraints (Eq. 1) to the ICRF2 defining sources. We used only sources with at least three 
sessions, ten delay observations, and a formal error ellipse major axis smaller than 5 nrad, and 
which were also in the ICRF2 defining list. 

Second, we did an after-the-fact verification that the NNR constraints indeed kept the positions 
from rotating with respect to the defining sources by estimating a 3-D rotation with respect to the 
ICRF2 defining sources. 

This second step was done using alternately diagonal-only parameter covariance and full covari- 
ance for the Celestial Reference Frame (CRF) being constrained. The reference, ICRF2, always 
used diagonal covariances because the off-diagonal correlations are not published. 

The results for three data sets are given in Table 1 with descriptions of the sets as follows: 

• The S/X data set included observations from October 1978 to September 2009 totalling 5.3M 
observations from 1709 sources. The NNR constraint used 267 ICRF2 defining sources. 

• The K-band data set [3] covers 2002 to 2009 with over 100,000 observations of 275 sources. 
The NNR constraint used 125 ICRF2 defining sources. 

• The X/Ka-band data set [2] covers 2005 to 2010 with over 10,000 observations of 387 sources. 
The NNR constraint used 153 ICRF2 defining sources. 


Table 1. 3-D Rotations (ga,s) estimated with vs. without correlations. 


CRF 

N Def 

Ri 

RRi 

R-2 

tr 2 

R:i 

°R 3 

S/X with corr. 

267 

1.4 

±5.4 

2.4 

±5.6 

2.0 

±4.4 

without corr. 

267 

-25.8 

±4.9 

9.4 

±4.9 

-1.6 

±4.3 

K with corr. 

125 

0.1 

± 6.8 

1.2 

± 7.4 

-0.1 

±5.2 

without corr. 

125 

-20.5 

±11.5 

-17.7 

±11.9 

-8.2 

±7.9 

X/Ka with corr. 

153 

0.3 

± 6.3 

2.9 

± 6.6 

- 0.1 

± 4.9 

without corr. 

153 

-35.1 

±17.9 

31.1 

±18.3 

38.6 

±11.2 


3. Discussion 

NNR Constraint Weighting: Given that the NNR constraint to the ICRF2 was unweighted, 
there may be problems caused by using a diagonal weighting for our after-the-fact check of rota- 
tional alignment. Our alignment check used the sum of the covariance of the frame under study 
and the diagonal ICRF2 covariance (note that the ICRF2 off-diagonal terms have not been made 
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available). For the three cases studied above, the X/Ka frame has the largest uncertainties, which 
are expected to dominate the combined X/Ka + ICRF2 covariance. The K-band uncertainties are 
somewhat larger than the ICRF2 and will partially dominate the combined covariance. However, 
because the S/X frame comes from almost the same data as the ICRF2, its uncertainties are com- 
parable to those of the un-inflated ICRF2. Once the ICRF2 diagonal covariances are inflated by 
a scale factor of 1.5 and then have a 40 //as noise floor RSS’ed into each position coordinate, the 
ICRF2 covariance dominates the diagonal uncertainties. 

These remaining nuances in weighting may explain why — even with full covariance — our after- 
the-fact check shows a few // as residual non-alignment. We have not yet studied these issues closely 
and can only say that the subject deserves more careful investigation. 

Nutation constraints: One important lesson learned in the course of our comparisons is that 
other constraints can distort the desired effect of the NNR constraint. In the case of our K-band 
solution, our a priori nutation model was the IAU 2000 standard of Matthews, Herring, Buffet 
(MHB) [6]. Because there was evidence external to our study that the MHB model had 100-150 
//as errors, we at first estimated nutation corrections with “weak” 200 //as a priori uncertainties 
on the nutation angles for each of twelve day-long sessions. Noting that a priori constraints act as 
pseudo-measurements, this effectively added 12 observations of the nutation pole each with 200 //as 
uncertainty. Thus, in retrospect, we realized that we had added a global constraint on the nutation 
pole of 200 //as / \/l2 = 58 //as which was strong enough to bias the 3-D rotational alignment of 
our celestial reference frame. This example is shared as a reminder that celestial frame alignment 
at the level of 10 //as requires careful attention to many details. 

4. Conclusions 

The celestial frame No-Net-Rotation constraint, soi x As* = 0, alters the a-5 full covariance 
produced by a CRF solution. The correlations contain significant information which is needed to 
correctly estimate 3-D rotations and associated sigmas. At the 5-10 //as levels of axial alignment 
precision found in current CRF work, these correlations must be accounted for in order to avoid 
statistically significant misalignment of a given CRF. 
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Figure 1. K-band (24 GHz) 12-session distri- 
bution of a-a correlations vs. arclength for all 
combinations of two sources. The correlations 
are systematically positive for short arcs and 
become negative for the longest arcs. These 
trends alter both the estimated rotations and 
their uncertainties. Color coding of 22.5° dec. 
bands: low to high (orange, red, green, blue, 
purple, black). 
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Figure 2. X/Ka (8.4/32 GHz) CRF’s distribu- 
tion of 5-d correlations vs. arclength. The cor- 
relations remain systematically positive over 
the full range of arc lengths, thus altering both 
the estimated rotations and their uncertain- 
ties. Color coding is the same as in the pre- 
ceding figure. Note the tendency for low dec- 
lination to be more strongly correlated. 


6. Appendix: Derivation of NNR Constraint Equation 

Given two sets of sources forming two frames, one frame is a priori assigned the role of the 
absolute reference and the other frame we wish to estimate from our data. Next, we establish 
the mathematical formula which describes the angular rotational offset between these two frames. 
Then we estimate the positions from data by least squares such that the rotational offset is zero. 
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Let soi be the reference position vector of the i th source, and let Sj be the estimated position. 
Denote the difference as As* = sot — s i ■ Suppose an approximately infinitesimal rotation, R(e ), is 
applied to s*. Then the rotated position is 

R(e) Si = Si + ex Si (7) 

We want the rotation R(e) that minimizes the mean square angular position difference: 

N src 

(A Q) 2 = —— Y (As; ~ e x s i) ' (As; - exs t ) ( 8 ) 

iV src ;= i 

We minimize with respect to variations in e.j by setting the derivative, 8 (A O) 2 / dej = 0 

N src Q 

Y lb - (As; - e x Sj)] • (Asi - e x Si) + (As* - e x Si) ■ [^7 (As; - e x st)] = 0 (9) 

;= l 3 

N src Q 

Y x Si) ■ (Asi - ex Si) + (As* - e x Si) ■ (-( £-e) x Si) = 0 ( 10 ) 

i= 1 ■? 

Since dot products commute, this can be re-written: 

-A svc 

Y (-2)(As l -exs l )-((i j e)xs i ) = 0 ( 11 ) 

i= 1 

We want the above derivative to be zero, d(AQ) 2 / dej = 0 when evaluated at e = 0, thus 

N src 

Y (A Si ) • ((£-e) x Si) = 0 for j = 1 , 2 , 3 ( 12 ) 

i= 1 

Using the general vector property v\ ■ (v 2 x v%) = V 2 ■ (v^ x v\) we have: 

N src 

Y • ( s i x As*) = 0 for j = 1 , 2 , 3 (13) 

i= 1 

Since is just the elementary unit vector on the j th axis, this represents three equations which 
can be combined into the vector equation: 

N src 

Y s i x A Si = 0 (14) 

i = 1 

Noting that s; x A Si = ( so; — A Si) x Asi = so; x Asi — A Si x A Si = so; x A.s* since for any vector, 
v, v x v = 0 , we have: 

N src 

Y So ; x Asi = 0 (15) 

i= 1 
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